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1. Introduction 

Physics analyses at the LHC will benefit if accurate predictions for background and signal 
processes become available. Arriving at such predictions often requires next-to-leading or- 
der (NLO) QCD computations. This is particularly true for multi-particle processes where 
the tree-level scattering amplitudes involve the strong coupling constant at a high power. 
In those cases, changing the renormalization scale often leads to 0(1) changes in the cross- 
section and more accurate predictions can only be obtained with NLO computations Q. 

The need for NLO corrections to processes with a vector boson and jets is particularly 
pressing. Corrections to vector boson + 1 jet processes and vector boson + 2 jet processes 
have been presented in refs. ||, [3|, |3], |j| and have been successfully compared with data in 
refs. H |7]]. The processes PP — > W/Z + N jets for N > 3 have a special importance. They 
constitute the principal background to a number of processes, such as top-pair production 
and i-channel single top production. In addition, PP — > W/Z + N jet production is an 
important source of jets + missing energy events, which is often regarded as a key channel 
in the search for physics beyond the Standard Model. 
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Techniques for NLO computations in the Standard Model in general and in QCD in 
particular are well developed. Traditional methods for NLO calculations are based on 
the observation that each Feynman diagram can be represented as a linear combination 
of tensor integrals. These tensor integrals can be reduced to scalar four-, three-, two- 
and one-point functions by exploiting Lorentz invariance; this procedure is known as the 
Passarino-Veltman reduction technique ||. 

While recent refinements of this procedure || 10, 11, |l3|, [l4|, |1| have transformed 



it into a powerful computational tool, there are two problems inherent in it. First, the 
number of diagrams in a particular process grows faster than N\ where N is the number 
of external particles. While processes with five or more external particles are rare at 
the Tevatron, the increase in energy and luminosity of the LHC makes consideration of 
processes with N > 5 particles phenomenologically mandatory. Second, in the course of 
the Passarino-Veltman reduction procedure for high-point functions, there are numerical 
instabilities related to the appearance of Gram determinants. The severity of this problem 
also increases with the number of external particles, because of the concomitant increase 
in the rank of the integrals. These two problems make the application of the Passarino- 
Veltman reduction technique to processes with more than five external particles highly 
non-trivial ju], |l7|, |l8|]. For example, currently there is not a single full process with six 
external particles for which NLO QCD corrections are known. 

While it may happen that traditional methods of one loop computations are able 
to overcome these problems |1S, 13], it is important to develop alternative solutions. One 
promising approach is the method of generalized unitarity that has been developed by Bern, 
Dixon, Dunbar and Kosower |2(J. Advances by Britto, Cachazo, Feng |pl|, [2^| allowed the 
development of analytic methods for the calculation of the full amplitude, including the 
rational part, using recursion relations [23, 24, 25]. A further recent advance by Ossola, 
Pittau and Papadopoulos p6| energized attempts to develop numerical procedures based 
on unitarity p7|, ||, ||]. 

A new computational scheme based on D-dimensional unitarity has been developed 
in Ref. |30|] . We will refer to this method as generalized D-dimensional unitarity. In 
Refs. [31, 32 1 the generalized D-dimensional unitarity method was further developed and 
was shown to be quite robust. In particular, it was explicitly demonstrated |3l]] that 
generalized D-dimensional unitarity is an algorithm of polynomial complexity where the 
evaluation time for one-loop pure gluonic amplitudes with iV external particles scales like 
N 9 . Moreover, it was also shown that generalized /^-dimensional unitarity can be applied 
to processes with massive fermions [32|. The results of these studies strongly suggest that 
generalized D-dimensional unitarity is an efficient computational algorithm which is now 
in a position to have a phenomenological impact. 

The goal of this paper is to make the first steps towards the application of generalized 
D-dimensional unitarity to phenomenology. We focus on the computation of virtual one- 
loop corrections to one of the important background processes at the Tevatron and LHC 
for which the one-loop corrections are still unknown - the production of the W boson in 
association with three jets. To this end, we have to consider one-loop corrections to the 
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processes 



0^u+d+g+g+g+ W + , 

0^u + d + Q + Q + g + W + , (1.1) 

and may assume, without loss of generality, that the quark Q does not couple to the 
W boson. We demonstrate that straightforward application of generalized D-dimensional 
unitarity allows us to compute all matrix elements required for the description of these 
complicated processes 1 . 

The paper is organized as follows. In the next Section, we summarize the salient 
features of generalized D-dimensional unitarity. In Section 3 we discuss the Dirac algebra 
and the choices of the polarization vectors in four- and higher-dimensional space-times. In 
Section 4 we describe the computation of all primitive amplitudes relevant for the process 
— > udgggW + . In Section 5 we focus on the amplitudes with four quarks, a gluon and a 
W boson. We conclude in Section 6. Numerical results for a specific phase-space point are 
collected in Appendix [A|. 



2. The method 



The method of calculation that we employ in this article is generalized /^-dimensional 
unitarity. The method relies on the observation []30| that one-loop scattering amplitudes 
in QCD can be fully reconstructed once tree-level scattering amplitudes are known for 
complex on-shell momenta of external particles, in, say, six- and eight-dimensions. The ne- 
cessity of knowing tree-level scattering amplitudes in higher-dimensional space-times stems 
from the fact that in QCD one-loop amplitudes are divergent and require regularization. 
Such regularization is conveniently done by continuing the dimensionality of space-time 
from four to 4 — 2e. At the end of the calculation, the limit e — > is taken, but ves- 
tiges of the regularization survive as particular finite contributions (rational terms) in the 
scattering amplitudes. It was pointed out in |}(J that the rational part of the amplitude 
can be determined by exploiting the dependence of residues of one-loop amplitudes on the 
dimensionality of space-time. Since this dependence is linear, it is sufficient to know these 
residues in two different space-time dimensions to reconstruct the residue as a function of 
D. 

For technical reasons, it is convenient to deal with scattering amplitudes where exter- 
nal particles are ordered and no permutations are allowed. Such ordering, for example, 
automatically fixes the flavors of all internal lines in the highest-level Appoint function 
that contributes to a particular A-particle ordered amplitude. It is well-known that such 



ordering can be achieved without sacrificing gauge-invar iance [ p3| , 34, [3Ej. For tree-level 
amplitudes, ordering of external particles appears naturally in color-ordered amplitudes. 
For one-loop amplitudes color ordering does not automatically lead to a complete ordering 
of all particles in the amplitude. To achieve this, color-ordered amplitudes are further 



In this article, we do not consider loop corrections with massive top quarks. Those contributions can 
be obtained along the lines described in E}^]. 
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decomposed into primitive amplitudes p^]. Those primitive amplitudes can be computed 



with the help of the color-stripped Feynman rules [36, 37, 38]. Note, however, that for a 
given primitive amplitude only color-charged particles are ordered while color-neutral par- 
ticles must be inserted in all possible locations to achieve a gauge-invariant result. For our 
purposes, this implies that the ordering of the W boson is not fixed and we have to account 
for all possible insertions of the W bosons between u and d quarks in a given primitive 
amplitude. 

To define a primitive one-loop amplitude, we employ the following set of rules: 

• we order all external particles that carry SU(3) color charge; 

• we draw a parent diagram with the direction of all fermion lines fixed such that 
the loop is always on the right-side of an upwards oriented fermion line. The order 
of the external particles is defined by reading the diagram clockwise. This defines 
left-handed primitive amplitudes 2 ; 

• for an iV-point scattering process, in general, the parent must be given by an one- 
particle irreducible iV-point function, represented by a diagram with N propagators 
in the loop. For some orderings it may happen that such a parent does not exist, in 
this case we draw the diagram by adding dummy lines; 

• we construct all possible cuts of a parent and we throw away all cuts that contain 
any dummy line; 

• we process each cut as required by generalized D-dimensional unitarity; tree-level 
on-shell amplitudes, needed for the computation of residues, are calculated using 
color-stripped Feynman rules. 

The parent primitive diagrams that are required for the calculation of W + 3 jet amplitudes 
will be presented later in the paper. 

The calculation of residues of primitive amplitudes requires the knowledge of tree- 
level amplitudes in six- and eight-dimensional space-time for complex momenta of external 
particles. The necessary matrix elements are constructed by employing Berends-Giele re- 
currence relations (3j|. Recall, that these recurrence relations connect off-shell currents 
of different multiplicities and with different particle content. The on-shell scattering am- 
plitudes are obtained from the on-shell limits of those currents. For the purposes of this 
paper, we need to employ currents with up to six on-shell external particles; a particular 
example is a fermionic current with three different fermion flavors and a gluon that con- 
tributes to some cuts of the u + d + W + + Q + Q + g scattering amplitude. We point 
out that, in a numerical program, it is possible to define those currents in a recursive way 
treating the number of external gluons as a parameter; currently, this is a necessary, (but 
not sufficient) prerequisite for the construction of fully automated computer codes for NLO 
QCD computations. 



2 There are also right-handed primitive amplitudes, where the loop is to the left of an upwards oriented 
fermion line. Since left and right primitive amplitudes are related, in this article we only present left 
primitive amplitudes |36|. 
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3. Dirac algebra, spinors and polarization vectors for gauge bosons 



3.1 Four-dimensional case 



For Dirac matrices it is convenient to use the Weyl representation where the 7-matrices 
are given by 
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Consider a massless fermion with momentum p = (E,p x ,p y ,p z ) and let p+ 
Solving the Dirac equation for massless quarks, we find the following solutions 



u\=i(p) 



( VpT \ 

(Px + iPy)/y/P+ 






u\=-\(p) 



( \ 



(Px - iPy)/y/P+ 

V ~VpT J 



(3.1) 
E + p z . 

(3.2) 



where A = ±1 refers to fermion helicity. Because p+ vanishes for E = —p z , the solution for 
the fermion moving in the — z direction requires care. Taking the limit, we arrive at 



u\=i(p) 



( \ 

s/2E 


V J 



u x =-i[p) 



( \ 


v/2~E 

V J 



(3.3) 



It is easy to see that in the massless case, the anti-particle solutions of the Dirac equation 
are related to the particle solutions so that v\(p) = u_\(p). 

The polarization vectors for massless gauge bosons in four dimensions are also well 
known. We present them here for completeness. For a gluon with momentum 



p = sin # cos </>, sin 6* sin^, cos 9), 



the polarization vector reads 

e x(p) = ~~/= (0; cos ® cos 4 
v2 



sgn(E)Xi sin 0, cos 9 sin <fi + sgn(E') Xi cos 0, — sin 9) . (3.4) 



In this paper we consider outgoing gluons; for this reason, all scattering amplitudes are 
computed with the complex conjugate vector e^(p). 

3.2 L>-dimensional case 

Generalized D-dimensional unitarity requires the knowledge of tree-level scattering am- 
plitudes in higher-dimensional space-time. To compute those amplitudes, we need D- 
dimensional polarization vectors for gluons, as well as spinors for fermions in D dimensions. 
Polarization vectors for gluons were discussed in detail in 31] and we do not repeat that 



discussion here. Weyl fermion spinors in higher-dimensional space-time are constructed as 
follows. 
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We construct a spinor solution for a fermion with light-like momentum p by using an 
auxiliary light-like vector n such that n ■ p 7^ 



Uj(jp,n) 



P 



V2p~ 



n 



Xj ( n ) > u j (P> ") = X; («) 



n 



(3.5) 



Here p = p^T 11 , where the summed index \i runs over D components, (the first four of which 
are the 0, x, y, z) and are the Dirac matrices in D dimensions. The index j specifies the 
spinor polarization states. We choose the D-dimensional, p-independent spinors Xj D \ n ) 
in such a way that 



2( D / 2 -i) 

3=1 



Xi 



n ) ® x\ D ^ ( n ) 



n. 



(3.6) 



In this case, it is easy to see that the Uj(p,n) spinors satisfy both the Dirac equation for 
massless fermions and the completeness relation 



2(0/2-1) 
^2 Uj(p,n) ®Uj(p,n) 
i=i 



pnp 
2p ■ n 



p. 



(3.7) 



We conclude that Uj(p,n) is a valid choice for on-shell fermion states. 

The above construction involves an auxiliary vector n and, for this reason is quite 
flexible. Having such a flexibility turns out to be important, especially since we have to 
construct on-shell spinors for complex momenta. We give a few examples below. 

We consider a D-dimensional vector n = (no,n x ,n y ,n z ,{n ie m_4\}), choose no = 
1/2, n z = 1/2 and set all other components to zero. Then, we need to find the spinors \ 
such that 



2(0/2-1) 



^2 Xj \n) ®x) '(») 



-A D )t 



n 



\ (To - r.) 



(3.8) 



Since ITo x y z are an block-diagonal [39], with "blocks" being 4x4 matrices, a L>-dimensional 
spinor is constructed by simple iteration of the four-dimensional construction. The four- 
dimensional spinors are given by 



(4) 

Xi 



/l\ 




W 



(4) 

x 2 



( \ 




V- 1 / 



(3.9) 



In six dimensions the eight-component spinors are choosen to be 



(6) 

Xi 



(6) 

x 2 




(6) 
X3 




(6) 
X4 




(3.10) 



The case D = 8 is a simple generalization of the above construction. 



- 6 - 



We now present two alternative procedures to define fermionic spinors which we employ 
when the particular choice of the vector n leads to numerical instabilities. This occurs for 
the on-shell momentum p = (po,0, 0,Po) since (p ■ n) = 0. To handle this case, we change 
the vector n to n = (1/2,0,0,-1/2,0/5-4) in the above formulas. However, even this can 
be insufficient. Indeed, note that a complex momentum p = (0,p x ,p y ,0) can be light-like. 
In this case, we need to choose yet another n. We can take n = (1, 1, 0, 0, Od-a) and choose 
the following four-dimensional spinors 

/ \ 

(4) 

x 2 



(4) 

Xi 



/1\ 

1 



w 



i 



(3.11) 



V- 1 / 



The higher-dimensional spinors are obtained from these four-dimensional solutions along 
the lines discussed above (see Eq. fl3.10D ). 



4. Processes with two quarks, a W boson and gluons 

In this section we consider the one- loop scattering amplitudes 0^u + d + (n — 2)g + W + . 
We refer to u as q and d as q and suppress the label of the W and its decay products 
in scattering amplitudes. We note that for a given primitive amplitude the W boson is 
inserted in all possible places when the diagram is traversed in a clockwise direction from 
q to q. 

4.1 Color decomposition of the amplitude 

At tree-level, the 0— > q + q + {n — 2) gluons + W scattering amplitude can be written as 



•A n (1<?) 2g, 3p, . . . , Tig) 



g n-2 (TM3)...r^) i ; i ir(l ? - 1 2 5 ;a(3) 5 ,.., ( T(n) 9 ) ) 



(4.1) 

where 5 n _2 is the permutation group of (n — 2) elements and A^ ee (lg, 2 q ;a(3) g , . . . , cr(n) g ) 
are color-ordered amplitudes. For all the amplitudes computed in this paper, we take 
the Wud interaction vertex to be —1^(1 — 7s)/2, so that neither electroweak couplings 
nor the Cabibbo-Kobayashi-Maskawa matrix elements are included. The W + decays to 
v(qi) + e + {q2); to account for this, we replace the polarization vector of the outgoing W 
by 



/ u u(g 1 )j^ ± v(q 2 ) 
1 J {qi+q2? ' 



7± 



(1±T5). 



(4.2) 



The choice of the polarization vector e_ corresponds to the W boson interactions in the 
Standard Model. The generators of the SU(3) color group are normalized as Tr(T a T b ) = 5 ab 
and satisfy the commutation relation 



[T a ,T b 



This normalization allows us to employ the color-stripped Feynman rules [36, 37 
calculate color-ordered scattering amplitudes. 



(4.3) 
58| to 
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Figure 1: Parent diagrams for primitive amplitudes A„ (lg, 3 g , ...m g , 2 q , m + l g ...n g ) and 
^j£,[i/2] ^ o g , ...m g , 2 q , (m + l) 9 ...n s ). All other parent diagrams that contribute to this primi- 
tive are obtained by considering all possible insertions of the W boson without changing relative 
ordering of quarks and gluons. The shaded circle stands for dummy lines. 



At one-loop, the color decomposition becomes more complicated. Using the color basis 
of Ref. [40 1, the one- loop scattering amplitude can be written as a linear combination of 
left primitive amplitudes 



yil-loop/-, 2 S n ) — n 

-^n I -*-<?> ^qi °gi ■ ■ ■ i ,L g I — y 



-p=2 a£S n - 2 

x (-1)M£(1 5 -, a(p) a{3) 2„ a(n) q , ...,a(p+ 1) ) (4.4) 



1 " * i i . c i i / ; i oi.ii '..,o \ 1 1. \ i ) i / > ■ — i . 

n-l 

i=lcrG5 n _2/S' n; j 



where for p = 2 the factor (T ■ ■ ■ T)i 2 n — > (T X2 T xl )i 2 n and forp = n the factor (F ■ ■ ■ F) X1X2 
5 X1X2 . In the second term S n -j = 1>j-\ is the subgroup of S n -2 that leaves Gr^ invariant. 
The color factors read 



Grg } (3,...,n) = N c (T as ■ ■ ■ T an 
Grg(3;4,...,n) = 0, 



12 



Grg } (3, . . . , j + 1; j + 2, . . . , n) = tr(T a * ■ ■ ■ T aj+1 ) (T^+ 2 • • • T a «)£ , j = 3, . . . , n - 2, 

Mi) 



GrJJ^S, . . . , n) = tr(T^ . . . T^) 6£ . (4.5) 

Parent diagrams for primitive amplitudes that involve two quarks and gluons are shown in 
Fig. | 

4.2 Numerical results for — ► qqgggW 

We have extended the Fortran90 program Rocket to include the computation of primi- 
tive amplitudes with quarks, gluons and gauge vector bosons. Rocket computes primitive 
amplitudes in the four-dimensional helicity scheme [^1] , 42 1. By default, the computation is 
done with double precision and, if a particular phase space point is deemed numerically un- 
stable, it is recomputed with quadruple precision using the package developed in Ref. J|3|] . 
The scalar integrals are evaluated using the QCDLoop package of ref. [H. 



-8- 



Note that, since Fortran90 supports recursive functions, implementation of Berends- 
Giele recursion for an arbitrary number of gluons is straightforward and indeed has been 
done in our program. Therefore, at least in principle, we can compute one-loop amplitudes 
for the process udW + + n gluons where n is an arbitrary number, We have checked that 
our program produces gauge-invariant one-loop amplitudes and correct 1/e 2 and 1/e poles 
for n up to ten. For n < 3, we have checked that our results agree with the known one-loop 



amplitudes for W + 1 and W + 2 jets H, 36, 45, B%. For the description of W + 3 jets, we 



require one-loop amplitudes with five external partons and this is what we focus on in the 
remainder of this paper. 

We present numerical results for the primitive amplitudes with seven external particles 
at a particular kinematic point considered in Ref. The momenta are chosen to be 

Pi(q) = — (1, cos a cos /3, cos a sin /3, sin a) , 

P2(q) = — (— 1, sin0, cos0sin0, cos0cos0) , 

Ps(g) = — (— 1, — sin0, — cos 0sin<^>, — cos 0cos <j)) , 

pM = i,o,o), 

Ps(g) = -(l,cos/3,sin/3,0) , 

P6(e + ) = ^(1, cos 7 cos /?, cos 7 sin /3, sin 7) , 

P7( v ) = ~Pl ~P2 ~P3 -PA ~P5 -Pe , (4.6) 
where [i = 1 GeV and 

rt 7T , 7T 7T 2-7T _ 37 , A 

= —, 4>=—, « = -, 7 = — , cos/3 = . (4.7) 

4 ' r 6 3 ' 3 ' 128 y J 

The momenta j?6 an d P7 are used to define the polarization vector of the W boson, eq.(ff .2j). 

Our results for unrenormalized primitive amplitudes A% and An'^ 1 ^ 2 ^ are summarized in 
the Appendix, see Tables |l|-|| and Tables |5||| respectively. We have checked that all primi- 
tive amplitudes have correct divergences and are gauge invariant. Moreover, we have tested 
the validity of our results for primitive amplitudes by reproducing a diagrammatic compu- 
tation of color-ordered amplitudes by taking appropriate linear combinations of primitive 
amplitudes. Finally, our program reproduces the results for the leading-color primitive 
amplitude A§ (lg, 2 q , 3 S , 4 S , b g ) computed recently in Ref. [46]. 



Next, we address the issue of the numerical stability of the computation. As was done 
earlier for similar studies of gluon amplitudes, we take care of numerical instabilities by 
performing computations with higher precision. Since higher precision slows the compu- 
tation, it is desirable to use it only for the phase-space points that suffer from numerical 
instabilities. The question we have to address therefore is how to detect numerical insta- 
bilities. To study this, we generate 10 5 random phase-space points using Rambo [47] with 



minimal constraints E± > lO -2 -^, |ry| < 3 and Ai? = y At/ 2 + A<p 2 > 0.4 and calculate 
the primitive amplitudes for — > qqgggW with double and quadruple precision. For each 
phase-space point, we can check whether or not the double precision computation of a 
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Figure 2: Accuracy for 2", 3", 4+, 5~) (left panel) and for A%(lf, 3",4+, 5 fl 

panel) for 10 5 randomly generated phase-space points. The raw double precision data as well as 
the result of numerical improvements are shown (see text for details). The inset shows the same 
plots in a linear scale. 



primitive amplitude reproduces the analytically known results for double and single poles 
in e and if the system of equations is solved with sufficient accuracy for each residue. To 
explain the latter test, we remind the reader that each residue is completely parameterized 
by a certain number of coefficients. We can check how well these coefficients are computed 
by choosing a random loop momentum, calculating the residue and checking how well this 
residue is obtained from the previously computed coefficients. We assign a relative error 
to each coefficient following the mismatch in this reconstruction. These errors are used to 
estimate the total error in the calculation of the primitive amplitude. By requiring that 
the relative precision in the poles and in the amplitude is better than 10~ 3 we find that 
around 0.3% of the points are recomputed in quadruple precision. 

After unstable points are recomputed with quadruple precision, the numerical insta- 
bilities are under control. This is demonstrated in Fig. ^ for two primitive amplitudes 
A?- (l|j~, 2~, 3~, 4+, 5~) and A§ (1^, 3~, 4+, 5~, 2~) where we show the number of events as 
a function of the relative accuracy eo defined as the absolute value of the difference be- 
tween double and quadruple precision results, divided by the quadruple precision result. We 
note, however, that the numerical stability of the amplitudes illustrated in Fig. ^ is generic, 
largely independent of the choice of the primitive amplitude and helicities of quarks and 
gluons. In fact, the two amplitudes considered in Fig. [2] are on the two sides of the spectrum. 
The leading-color amplitude A^ (I^~, 2~, 3~,4+, 5~) has the minimal number of cuts, since 
the W boson can only be inserted in one place, between lq and 2 q . On the contrary, the 
amplitude A\ (1^, 3~, 4+, 5~, 2~) has the maximal number of cuts since the W boson can 
be inserted in four different places. Thus, among all primitive amplitudes with two quarks 
and three gluons, maximal computational effort is required for A^ (ljj", 3~, 4+, 5~, 2~) so 
that the issues of numerical stability may be expected to be worst in this case. 

We expect that further optimization of the procedure for identifying unstable points 
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may be required to arrive at an optimal compromise between numerical accuracy and speed 
of the code. For instance, with an arbitrary precision package such as that of ref. [p3|j . 
one can design a procedure where instead of using fixed quadruple precision for unstable 
points, the number of digits in the higher precision calculation is established according to 
how unstable the point is. We plan to study these issues in the future. 

Finally, we remark on the CPU-time required for the evaluation of one-loop amplitudes. 
We find that it takes about 45-50 msec to evaluate the leading-color primitive amplitude 
(lg, 2q, 3 9 , 4 g , 5 g ) on a computer with 2.33 GHz Pentium Xeon processor using the intel 
fortran compiler. This is comparable to the evaluation times for six-gluon primitive am- 
plitudes |HJ]. For more complicated primitive amplitudes, the number of cuts increases 
and the evaluation times scale accordingly. For example, for A§ (lq, 3 g , 4 ff , 5 g , 2 q ), the case 
with the maximal number of cuts, it takes about 160 msec to compute a single primitive 
amplitude. 



5. Processes with two quark pairs, a W boson and a gluon 
5.1 Color decomposition of the amplitude 

We now turn to processes with two quark pairs, the W boson and a gluon, — > u + d + 
Q + Q + W + g. We will again refer to u as q and to d as q. To construct color-ordered 
primitive amplitudes, we assume that the W boson can not couple to the quark Q. The 
color decomposition of the four-quark and one-gluon amplitude at tree-level reads 



(rna 5 \ ijr 13 Dtree , 1 (rna 5 \ lu iz Dtree 
l J )i A °i 2 ^7;1 >ii °U a 1\2 



(5.1) 



in an obvious notation where 05 is the color index of the produced gluon. The color 
decomposition of the four-quark and one-gluon amplitude at one loop reads 



»4 *2 

Each of these one-loop color-ordered amplitudes can be further written as 

i = 1,2,3,4, 



(5.2) 



(5.3) 



to separate the diagrams with a closed fermion loop from the other ones. The amplitudes 
Bj~\ and B^ 1 ^ can be written as linear combinations of primitive amplitudes. Those 
primitive amplitudes are shown in Figs. [3| and ||. 
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Q' 




Figure 3: Prototype parent diagram for primitive amplitudes with four quarks, a W boson and 
a gluon which contain a closed fermion loop. The gluon can be inserted in four possible ways into 
the prototype graph, leading to four primitive amplitudes. Note that the W only couples to q. The 
solid blobs denote the dummy lines introduced in section 2. 



For the amplitudes with an additional closed fermion loop we find 



R [l/2] 
D 7;\ 

R [l/2] 
7;2 

o[l/2] 
-°7;3 

o[l/2] 
-°7;4 



-4 1/2] (l g -,5 9 ,4 Q ,3Q,2 g ), 
-4 /2] (l g -,4 Q ,3Q,2 g ,5,), 
-A [ l /2] m,4 Q ,3 Q ,5 g ,2 q ), 
-4 /2] (l ? ,4 Q ,5 9 ,3Q,2 g ). 



(5.4) 



The three classes of primitive amplitudes that we need to consider for four-quark 
processes without closed fermion loops are shown in Fig. |], 



B 



[i] 

7;i 



(5.5) 



Amplitudes from each class are written as linear combinations of primitives amplitudes. 




Figure 4: Prototype parent diagrams for primitive amplitudes with four quarks, a W boson and 
a gluon for classes a, b and c. The gluon can be inserted in four possible ways into any of the 
prototype graphs, leading to four primitive amplitudes in each case. For the diagrams in class b, 
all possible insertions of a W boson to a given parent primitive should be considered; note that the 
W only couples to q. The solid blobs denote the dummy lines introduced in section 2. 
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For class a, we find 

4!i* a = (l - 4 ] ' a (l,-,2 9 ,3Q,4 Q ,5 a ) - ^ (-4 1] >° (1,-,5 s ,2 ? ,3q,4q) (5.6) 

-4 ha (Iff. 5 9) 2„ 4 Q , 3q) - 2„ 5 g , 3q, 4q) - Af a (l q , 2 g , 5 g , 4 Q ,3 Q ) 

-Af a (l,-, 2„ 3q, 5 S , 4q) - 4 1] ' a 2„ 4 Q , 5 fl , 3q) + Af a 2„ 4 Q , 3$, 5 g )) , 
4!i' a = ^ 11 ' a (^2 g ,5,,4 Q ,3Q)-4 1 ' a (l g -,2 g ,4 Q ,5 9 ,3Q)+4 1] ' a (l,-, 2„ 4 Q , 3q, 5 fl ) 

(4 ] ' a (iff, 5 9 , 2„ 3 , 4q) + 4' a (l,-, 5 g , 2„ 4q, 3q)) , (5.7) 

= (i - ^) 4 1I,B (1,-,2 9 ,5 s ,3q,4 ) - ^ (-4 1] ' a (1,-,5 s ,2 ? ,3q,4q) (5.8) 
-^(I^^^Sq)-^ 1 ^ 

-4 11 ' a (l ? -,2 ? ,3 ,4 Q ,5 fl )-4 11 ' a (l ? -,2 ? ,4 ,3Q,5 a )+4 11 ' a (l ? -,2 ? ,5 g ,4 ,3 )) ) 
4!l' a = -4 1,a (^5„2 9 ,4 Q ,3 ) - 4 ] ' a m,2 q ,5 g AQ,S Q ) ~Af a (1,-,2„4q,3q,5 p ) 

" W (1 " 2 " 4<3 ' 5ff ' ^ + (1 " 3<3 ' 59 ' 4Q) ) ' (5 ' 9) 

Note that in this formula there are amplitudes where fermions and anti-fermions alternate 
(qqQQ) and amplitudes where this is not the case (qqQQ)- The latter can be reduced to 
the former using the following C-parity relation 

/\l 1 h a f-t_ 9 0A3 \ — (_-\\n+i A 1 }' - f-\_ 9 A Xi ^ \ 

L \ L qi ■ ■ ■ > Z <?> • • • > 4 Q ' • • • ' 6 Q ' • • ■) — V L ) A L \ 1 Q^ ■ ■ ■ ' Z <?' • • • 4 Q ' • • • ' J Q ' • • -J • 

(5.10) 

Here n is the number of external gluons sandwiched between the Q and Q spinors. 
For classes b and c we obtain 

M!P = ^4 11,b (1* 5„ 4^30,2,), 

= -^5 K 11 ' 6 (Iff, 5 P , 4q, 3q, 2,) + Af b (1„ 4 Q , 5 fl , 3 , 2 ? ) + A™' 6 (l„ 4 Q , 3 Q , 5 P> 2,)) , 

fi S' 6 = ^24 1I,6 (lff.4«.3 <5 ,5 ff> 2,) > 

4! 6 = -^(iff.Sfl.^.Sg^,) -4 1] ' 6 (l ? ,4 Q ,3Q,5 fl ,2 9 ) -4 1] ' 6 (l 9 -,4 Q ,3Q,2 9 ,5 fl ) 

" ( X " ^2) 4 ] * fe (Iff, 4q, 5 ff , 3q, 2,) . (5.11) 

B S' C= i^4 1,C (lff,5 5 ,4 Q ,3 a ,2 ? ), 

B W- C = -4' c (1 5 ,5 9 ,4q,3q,2,) -4 1] ' c (l 9 -,4 Q ,5 9 ,3Q,2 g ) -4 ] ' c (1 5 ,4 q ,3q,5 9 ,2 9 ) 
"(l-^) 4 11,C (^4Q,3 ,2 g ,5 9 ), 
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V) 
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10 




-14 -12 -10 



Figure 5: Accuracy for ^4^'' a (l^, 2~, 5+, 3q, 4 J) amplitude for 10 5 randomly generated phase- 
space points. The raw double precision data as well as the result of numerical improvements are 
shown. The inset shows the same plots in a linear scale. 



-8 -6 -4 
l°9io( e o) 



B W,c 

-°7;3 



B, 



1 



7:1 



4 11,C (^4q,3 ,5 9 ,2 9 ), 



1 



[l],c 



(5.12) 



^(AW' c (l ? - ) 5 gj 4 Q ,3Q,2 ? )+^' c (l $ ,4Q ! 3 Q ,5 g5 2 g )+^' c (l $ ,4 Q ,3Q,2 ? ,5 



5.2 Numerical results for — > qqQQgW amplitudes 

We now present numerical results for — ► qqQQgW amplitudes. We use the kinematic 



point in Eq. (4.6), where the momentum is the momentum of Q and the momentum 
P4 refers to Q. The results of the calculation are given in Tables [7-11. The results for 
primitive amplitudes given in these Tables are sufficient to obtain numerical results for 
the color-ordered amplitudes B 7 -i, i = 1, ... 4 using the equations given in the previous 
subsection. In fact, we have checked our results for the primitive amplitudes by calculating 
the color-ordered amplitudes Bj-i, i = 1, ... ,4 diagrammatically and verifying that linear 
combinations of the primitive amplitudes computed with Rocket reproduce results of the 
diagrammatic computation. The numerical stability of a four-quark amplitude is illustrated 
in Fig. H It is very similar to the numerical stability of the amplitudes with two quarks 
and three gluons discussed earlier in detail and to other four-quark amplitudes. 



6. Conclusions 

In this paper we have shown that a straightforward application of the generalized D- 
dimensional unitarity method proposed in ref. |3(| allows us to compute all one-loop scat- 
tering amplitudes required to describe the production of a W boson in association with 
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three jets at hadron colliders. We observe satisfactory performance in terms of numerical 
stability and required run times. We feel confident that the results of this paper provide 
a solid foundation for computing the one-loop virtual corrections to the production of the 
W + 3 jets in hadron collisions. 

On a more general side, the current version of Rocket computes one-loop amplitudes for 
processes — > n gluons, — ► qq+n gluons, — > qqW+n gluons and — > qqQQW+1 gluon. 
It is straightforward to extend the program to include similar processes with the Z boson 
and processes with massive quarks — > it + n gluons. This list is a testimony to the power 
of the method and indicates that the development of automated programs for one-loop 
calculations may finally be within reach. 
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A. Numerical results 



In this Appendix we present the numerical results for various primitive amplitudes that 
we employ for the computation of W + 3 jet processes. Numerical results are presented for 
the phase-space point given in Eq.( |4.6| ) for various helicities of the external particles. For 
convenience, we present the results for the ratio of a one-loop primitive amplitude and the 
corresponding primitive tree-level amplitude defined as 3 

r % 2 3 4 5 6 7)- 1 feWl c, _ r(l + g )r(l-6) 2 

r L (1, 2, 3, 4, 5, b, 7) - ^ ^ ^ ^ ^ ^ ?) , c r - (47r)2 „ er(1 _ ^ , (A.l) 

where in this appendix we always indicate explicitly the dependence on the lepton momenta 
from the W decay. 

A.l Numerical results for — > qqgggW amplitudes 



Helicity 


1/6 2 


1/6 


e° 


At-(l+2-3+4+5+6±77-) 
4 1] (1^2-3+4+5+6±77-) 


-4.00000 


-10.439578 -i 9.424778 


-0.006873 + i 0.011728 
5.993700 - i 19.646278 


A tr -(l+2-3+4+5-6±7 r ) 
4 1] (^2-3+4+5-6±77-) 


-4.00000 


-10.439578 -i 9.424778 


0.010248 -% 0.007726 
-14.377555 -i 37.219716 


^(l±2-3-4+5+6±7D 
rf ] (l±2-3-4+5+6±7r-) 


-4.00000 


-10.439578 -i 9.424778 


0.495774 - il. 274796 
-1.039489 - i 30.210418 


^-(1+2-3-4+5-6^77-) 
*-2 ] (lf^3-4+5-^7T) 


-4.00000 


-10.439578 -i 9.424778 


-0.294256 -i 0.223277 
-1.444709 - i 26.101951 



Table 1: The primitive tree-level amplitude ^4 troo (lg, 2 q , 3 g , 4 9 , 5 S , 6^, 7;) and the ratio 
rj} (lq, 2 q , 3 g , Ag, 5 g , 6r, 7i) of a primitive one-loop amplitude to the primitive tree-level amplitudes 
for various helicities. 



3 Note that A tTce (l, 2, 3, 4, 5, 6, 7) denote here primitive tree-level amplitudes, as opposed to the color 
ordered amplitudes ^4j; ree (lq, 2 q ; a (3) a , . . . , cr(n) g ) appearing e.g. eq. ( |4.l| ), where the quarks are separated 
by a semicolon from the gluons. 
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Helicity 


1/6 2 


1/6 




^ tree (l±3+2-4+5+6± 7,-) 
»"P(l3"3+2-4+5+6±V) 


-3.00000 


-8.676830 - i 6.283185 


-0.005446 + i 0.009804 
-1.423339 - i 14.443863 


A*"-(l+3+2-4+5j6±7T) 

h 1 

^(^3+2-4+5-6^7^ 


-3.00000 


-8.676830 - i 6.283185 


0.000364 + i 0.004550 
-11.406265 - i 16.485295 


A trcc (l+ 3-2-4+5+6^7,-) 
4 1] (^3-2-4+5+6+7 r ) 


-3.00000 


-8.676830 - i 6.283185 


0.341643 -i 0.310960 
-5.430180 -i 21. 180247 


A^e ( i + 3-2-4+5-6+7 r ) 
rt 1] (l+3-2-4+5-6+7 r ) 


-3.00000 


-8.676830 -i 6.283185 


0.024966 -i 0.156703 
-4.868668 -i 21.036597 



Table 2: The primitive tree-level amplitude A trcc (lq, 3 S , 2 q , A g , 5 g , 6f, 7i) and the ratio 
r^'(lg, 3 9 , 2 q , A g , 5 g , 6,, 7;) of a primitive one-loop amplitude to tree-level primitive amplitude for 
various helicities. 



Helicity 


l/e 2 




e° 


^-(1+3+4+2-5+6+7,-) 
4 11 (1^3+4+2-5+6+7,-) 


-2.00000 


-7.835662 -i 3.141593 


-0.005563 - i 0.030746 
13.662096 -^25.637707 


A tree (l+3+4+2-5-6+7,-) 
^(^3+4+2-5-6+7,-) 


-2.00000 


-7.835662 -^3.141593 


0.022677 + i 0.085524 
-9.177581 -i 16.265480 


^-(1+3-4+2-5+6+7,-) 
^(^3-4+2-5+6+7,-) 


-2.00000 


-7.835662 -^3.141593 


-0.098988 + t 1.958409 
-12.140461 -i 15.924761 


^-(1+3-4+2-5-6+7,-) 
^(^3,-4+2-5-6+7,-) 


-2.00000 


-7.835662 -^3.141593 


-0.283565 +i 0.841833 
-13.465828 -i 13.730719 



Table 3: The primitive tree-level amplitude A trcc (l q , 3 9 , A g , 2 q , 5 S , 6j, 7i) and the ratio 
r^'(lg, 3 S , A g , 2 q , 5 g , 6f, 7/) of a primitive one-loop amplitude to the primitive tree-level amplitude 
for various helicities. 
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Helicity 


1/6 2 


1/6 


e° 


A*~(l+3+4j5+2-6±7T) 
^(^3+4+5+2-6^7,-) 


-1.00000 


-3.334232 - i 0.000000 


0.017883 + i 0.009214 
11.973924 -i 8.033958 


^-(1+3+4+5-2-6+7,-) 

r-i ] 

^(^3+4+5-2-6^7,-) 


-1.00000 


-3.334232 + i 0.000000 


-0.033289 - i 0.082348 
-1.783190 + i 1.552944 


A trcc (l 9 +3-4+5+2-6+ 7,-) 
4 11 (^3-4+5+2-6±7,-) 


-1.00000 


-3.334232 - i 0.000000 


-0.738428 - i 0.372652 
-5.654597 + z 0.276608 


^(1^3-4+5-2-6+7,-) 
4 11 (^3-4+5-2-6±7,-) 


-1.00000 


-3.334232 +i 0.000000 


0.552856 -i 0.461853 
-6.461431 +i 1.451815 



Table 4: The primitive tree-level amplitude A trcc (lg, 3 9 , A g , 5 S , 2 q , 6j, 7i) and the ratio 
r^'(lg, 3 9 , A g , 5 fl , 2 g , 6,, 7;) of a primitive one-loop amplitude to the primitive tree-level amplitude 
for various helicities. 



Helicity 


1/e 2 




e° 


A trcc (l 9 +2-3+4+5+6+ 7,-) 






-0.006873 + i 0.011728 


r[ 1/2l (l+2-3+4+5+6+7,-) 


0.00000 


0.000000 -t 0.000000 


-6.001512 - i 26.601839 


A^ c (l± 2" 3+4+5-6 r +7,-) 
r[ 1/2] (l+2-3+4+5 9 6+7,-) 


0.00000 


0.000000 -i 0.000000 


0.010248 -i 0.007726 
7.227836 -i 4.090839 


A trcc (l 9 +2-3-4+5+6f7,-) 






0.495774 - i 1.274796 


r[ 1/2l (l+2-3 9 -4+5+6+7,-) 


0.00000 


0.000000 +i 0.000000 


0.096288 -i 0.114398 


^ tree (l±2-3j4+5-6 r + 7,-) 
4 1/21 (l?2-3-4+5-6±7 r ) 






-0.294256 -^0.223277 


0.00000 


0.000000 +i 0.000000 


0.164410 - i 0.134601 



Table 5: The primitive tree-level amplitude A trec (lg, 2 q , 3 g , 4 S , 5 S , 6f, 7i) and the ratio 
r L^ 2 ' 2(j , 3 9 , A g , 5 g , &i, 7i) of a primitive one-loop amplitude to the primitive tree-level amplitude 
for various helicities. 



- 18 - 



Helicity 


1/e 2 






A trcc (l±3+2"4+5+6^ 7,-) 
r[ 1/2l (l+3+2-4+5+6±7 r ) 


0.00000 


0.000000 -i 0.000000 


-0.005446 + i 0.009804 
-0.127241 - i 1.316987 


A^(l+3+2-4+5j6±7f) 
r[ 1/21 (1^3+2-4+5-6±7r-) 


0.00000 


0.000000 +i 0.000000 


0.000364 + i 0.004550 
0.000000 + i 0.000000 


A tree (l±3-2-4+5+6f 7f) 
rg^(l+3j2-4+5+6f 7T) 


0.00000 


0.000000 +i 0.000000 


0.341643 -i 0.310960 
-0.020783 -i 0.001593 


A tree (l±3" 2-4+5" 6^ 7,-) 
r[ 1/2l (1^3-2-4+5 g -6±7 r ) 


0.00000 


0.000000 -i 0.000000 


0.024966 -i 0.156703 
0.000000 - i 0.000000 



Table 6: The primitive tree-level amplitude A tree (lq, 3 S , 2 q , A g , 5 S , 6j, 7i) and the ratio 
r L^ 2 ' (lg> 3 g , 2 g , 4 g , 5 g , 6j, 7i) of the primitive one-loop amplitude to the primitive tree-level am- 
plitude for various helicities. 



- 19 - 



A. 2 Numerical results for — > qqQQgW amplitudes 

Numerical results for amplitudes with two quark pairs, a W boson and a gluon are presented 
in Tables below. 



Helicity 


1/e 2 


1/e 

/ 


e° 


A*-(l+5+2-3g 4J6±7D 
rt 1] ' a (l±5+2-354+6±7r) 

J-' v ^ y h (erf "^z I 1 


-2.00000 


-1.906388 -i 6.283185 


-1.347977 -i 0.593626 
12.513239 - i 16.727811 


r[ 1] ' a (l±5-2-354+ 6±7f) 


-2.00000 


-1.906388 - i 6.283185 


-0.570749 - i 0.316836 
12.452841 - i 14.482266 


^e (1 + 2 - 5 +3- 4 + 6 + 7r) 
rt 1] ' a (l±2-5+354+6±7r) 


-3.00000 


-6.083408 - i 3.141593 


0.734834 - i 0.360895 
9.466663 - i 4.461718 


A*"-(l+2-5-3^4+^7T) 
rt 1] ' a (l±2-5j334+6±7f) 


-3.00000 


-6.083408 - i 3.141593 


0.421057 + i 0.463392 
7.337316 - i 3.094966 


^ r -(l+2-3^5+4+6±7-) 
4 11,a (ir2-3 5+4+6±7r) 


-2.00000 


-1.906388 - i 6.283185 


-0.288607 -i 0.043034 
14.470591 - i 16.520704 


^e (1 + 2 -3- 5 - 4 J 6 + 7r) 
rW' a (l±2 g "3 5-4+6±7 r ) 


-2.00000 


-1.906388 - i 6.283185 


-0.148841 + i 0.001847 
10.478262 - i 16.315474 


^- (1 + 2 - 3 - 4 j 5 + 6 + 7r) 


-3.00000 


-5.946351 - i 9.424778 


0.901749 + i 0.997556 
10.012255 - i 29.539947 


^ ree (l g +2-334+5 9 -6f7 r ) 
rW' a (l±2-33 4+5-6±7 r ) 


-3.00000 


-5.946351 -i 9.424778 


0.298533 -i 0.147741 
9.462997 - i 25.465941 



Table 7: The primitive tree- level four-quark amplitudes A trec and the ratio r L of the one-loop 
class a primitive amplitude to the corresponding primitive tree-level amplitude for various helicities. 
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Helicity 


l/e 2 




e° 




-2.00000 


-3.109446 + i 0.000000 


-1.347977 - i 0.593626 
7.650907 -i 1.630983 


^-(l±5-2-4+ 3~ Q 6±7i) 
4 1] ' a (^5-2-4j3 5 6+7 r ) 


-2.00000 


-3.109446 + z 0.000000 


-0.570749 -i 0.316836 
7.388565 +i 2.000057 


^ ree (l ? +2-5+4+3Q6±7r) 
rP- B (lj2-5+4+3 5 ^7r) 


-3.00000 


-6.730261 -^3.141593 


0.446227 -i 0.403929 
8.835307 -i 10.218599 


^ tree (l g +2-5-4t3 Q 6f 7-) 
rW' a (l±2-5-4j3Q6±7 r ) 


-3.00000 


-6.730261 -^3.141593 


0.272217 + i 0.464577 
2.417473 - i 12.686072 


^ ree (l g +2-4j5+3Q6±7 r ) 


-2.00000 


-3.109446 + i 0.000000 


0.288607 + i 0.043034 
11.636607 + i 2.670357 


^(l±2-4j5-3 5 6±7 r ) 

4 11,a (^ 2-43 5-3^-) 


-2.00000 


-3.109446 + i 0.000000 


0.148841 -i 0.001185 
3.098383 -i 0.972210 


A*"-(l+2-4+3g5+^V) 
rW' a (l±2-4j3 5 5+6±7 r ) 


-3.00000 


-6.502556 -^3.141593 


0.613143 + ^0.954522 
9.563123 - i 12.742650 


^(l±2-4j3Q5-6±7 r ) 
4 1] ' a (^2-4j3 Q 5-6±7 r ) 


-3.00000 


-6.502556 -^3.141593 


0.149692 -i 0.146556 
7.832425 -^8.905177 



Table 8: The primitive tree-level four-quark amplitudes ^4 tree and the ratio a of the primitive 
one-loop class a amplitude to the corresponding primitive tree-level amplitude for various hclicities. 
Note the different ordering of the quarks Q and Q compared to Table 7. 



- 21 - 



Helicity 


l/ e 2 


l/e 


e° 


^ tree (l^ 5+ 4+ 3"= 27 6^ 77) 
^'^5+ 4+ 37,27 6,+ 77) 


-1.00000 


-3.334232 + i 0.000000 


-0.901749 - i 0.997556 
-8.027321 - i 0.968722 


A tree (l+ 57 4+ 37, 2" 6+ 77) 


-1.00000 


-3.334232 +i 0.000000 


-0.298533 + i 0.147741 
-4.397007 -i 3.109084 


^(1+4+ 5+ 3^27 6+ 77) 
rP' 6 (lj4+5+3^2-6j!-V) 


-1.00000 


-3.334232 + i 0.000000 


0.288607 + i 0.043034 
-5.997531 - i 0.856042 


^-(1+4+ 57 3^27 6+ 77) 


-1.00000 


-3.334232 +i 0.000000 


0.148841 -i 0.001185 
-6.462353 - i 0.674696 


^-(1+4$ 3^5+27 6+ 77) 
^1^37,5+27 6+77) 


-1.00000 


-3.334232 - i 0.000000 


-0.734834 + i 0.360895 
-6.039221 - i 0.200173 


^(1+4+3^57 27 6+77) 
rE ] **W4+3^5-2-6f 7f) 


-1.00000 


-3.334232 +i 0.000000 


-0.421057 - i 0.463392 
-6.632788 -i 0.178166 


^-(1+37,4+27 5+6+77) 
3^4+27 5+ 6+77) 


-2.00000 


-7.835662 -i 3.141593 


1.347977 + t 0.593626 
-12.743438 - i 16.092839 


A*"-(1+3q4J2-5J^7T) 
^1^4+ 27 57 6+77) 


-2.00000 


-7.835662 -i 3.141593 


0.570749 + i 0.316836 
-14.463179 - i 14.071065 



Table 9: The primitive tree- level four-quark amplitudes A trcc and the ratio r^' b of the one-loop 
class b primitive amplitude to the corresponding primitive tree-level amplitude for various hclicities. 
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Helicity 


l/ e 2 


1/e 


e° 


,4 tree (l + 5+ 4+ 3. 27 6+ 77) 
y 1 9 Q Q 1 I I > 

^(^5+ 4+ 3^27 6+ 77) 


-1.00000 


-3.826559 -i 0.000000 


-0.901749 - i 0.997556 
-10.068901 - i 0.284089 


,4 tree (l+ 57 4+ 3"= 27 6+ 77) 
^'W^ 4+ 3^27 6+ 77) 


-1.00000 


-3.826559 -i 0.000000 


-0.298533 + i 0.147741 
-9.118256 -i 0.316181 


^-(1± 4+ 57 3^27 6+ 77) 
^'W 4$ 57 3^27 6+ 77) 


-2.00000 


-5.954375 - i 3.141593 


0.288607 + i 0.043034 
-6.343258 - i 9.293058 


^-(1±4J57 3327 6±77) 
^ 1] ' C (lf4j 57 3^2- 6±7 r ) 


-2.00000 


-5.954375 -i 3.141593 


0.148841 -i 0.001185 
-8.567533 -t 11.440611 


^-(1+4$ 3^5+27 6+ 77) 
^(13-4+3^5+2-^"^) 


-1.00000 


-3.826559 + i 0.000000 


-0.734834 + i 0.360895 
-9.495931 + a 0.100657 


4+ 3^57 27 6+ 77) 
^ 11 ' C ( 1 y ±4 Q 3 5 9 " 2 y 6 fV) 


-1.00000 


-3.826559 -i 0.000000 


-0.421057 - i 0.463392 
-9.538700 -^0.299195 


A*«(l+3g4j2-5+^V) 
4 11 ' C (ly ± 3Q4j27 5+6 + 7 r ) 


-1.00000 


-3.826559 + i 0.000000 


1.347977 + i 0.593626 
-9.696279 + i 0.000000 


A^(l± 37, 4± 27 5^6+77) 
4 11 ' C (ly ± 3Q4j27 5-6±7 r ) 


-1.00000 


-3.826559 + i 0.000000 


0.570749 + i 0.316836 
-9.696279 + i 0.000000 



Table 10: The primitive tree-level four-quark amplitudes A trcc and the ratio r L of the primitive 
one-loop class c amplitude to the corresponding primitive tree-level amplitude for various hclicitics. 
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Helicity 


l/ e 2 






^ tree (l^ 5+ 4+ 3^ 27 6^ 77) 

J 1 /2]/ 1 + r + 4 + o- o- «+ 7-^1 


0.00000 


-0.666667 - i 0.000000 


-0.901749 - i 0.997556 
-2.527058 - i 0.104842 


A tTee (l± 57 4+ 3"= 27 6± 77) 
4 1/2l ( 1 y ± 57 4+3^27 6±7 r ) 


0.00000 


-0.666667 + i 0.000000 


-0.298533 + i 0.147741 
-2.366982 +i 0.089047 


^(l±4+5+3^27 6±77) 
4+ 5+ 37,27 6,+ 77) 


0.00000 


-0.66666 + i 0.000000 


0.288607 + i 0.043034 
-2.656572 + i 0.000000 


^-(1± 4$ 57 3^27 6+ 77) 
4+57 3^ 27 6+ 77) 


0.00000 


-0.666667 + i 0.000000 


0.148841 -i 0.001185 
-2.656572 +i 0.000000 


^-(1+4$ 3^5+ 27 6+77) 
4 1/2] (44+3^5+27 6±7 z -) 


0.00000 


-0.666667 + i 0.000000 


-0.734834 + i 0.360895 
-2.887087 - i 0.164880 


A^(l±4j3^57 27 6±77) 
4 1/2l (l^S3^57 27 6±7 r ) 


0.00000 


-0.666667 + i 0.000000 


-0.421057 - i 0.463392 
-2.749917 + i 0.137009 


^e (1 + 4 +3- 2 - 5 + 6 ± 7r) 
(1^37, 27 5+ 6+77) 


0.00000 


-0.666667 + i 0.000000 


1.347977 + i 0.593626 
-2.662151 +i 0.000000 


^-(l±4+3^27 57 6±7 r ) 
4 1/2l (l?4S3 5 27 57 6±7 r ) 


0.00000 


-0.666667 + i 0.000000 


0.570749 + i 0.316836 
-2.662151 +i 0.000000 



Table 11: The primitive tree- level four-fcrmion amplitude ^4 tree and the ratio r L of the primitive 
one-loop amplitude to primitive tree-level amplitude for various helicities. 



- 24 - 



References 

[1] Z. Bern et al. [NLO Multileg Working Group], arXiv:0803.0494 [hep-ph]. 

[2] R. K. Ellis, G. Martinclli and R. Pctronzio, Nucl. Phys. B 211, 106 (1983). 

[3] P. B. Arnold and M. H. Reno, Nucl. Phys. B 319, 37 (1989) [Erratum-ibid. B 330, 284 
(1990)]. 

[4] W. T. Giele, E. W. N. Glover and D. A. Kosowcr, Nucl. Phys. B 403 (1993) 633 
[arXiv:hcp-ph/9302225]. 

[5] J. Campbell and R. K. Ellis, Phys. Rev. D 65, 113007 (2002) [arXiv:hep-ph/0202176]. 

[6] T. Aaltoncn et al. [CDF Collaboration], Phys. Rev. D 77, 011108 (2008) [arXiv:071 1.4044 
[hep-ex]]. 

[7] T. Aaltoncn et al. [CDF - Run II Collaboration], Phys. Rev. Lett. 100, 102001 (2008) 
[arXiv:071 1.3717 [hep-ex]]. 

[8] G. Passarino and M. J. G. Veltman, Nucl. Phys. B 160, 151 (1979). 

[9] A. I. Davydychev, Phys. Lett. B 263 (1991) 107. 

[10] G. Duplancic and B. Nizic, Eur. Phys. J. C 35 (2004) 105 [arXiv:hcp-ph/0303184]. 

[11] W. T. Giele and E. W. N. Glover, JHEP 0404 (2004) 029 [arXiv:hcp-ph/0402152]. 

[12] R. K. Ellis, W. T. Giele and G. Zanderighi, Phys. Rev. D 73 (2006) 014027 
[arXiv:hcp-ph/0508308]. 

[13] A. Dcnncr and S. Dittmaicr, Nucl. Phys. B 734, 62 (2006) [arXiv:hcp-ph/0509141]. 

[14] A. van Hamcrcn, J. Vollinga and S. Weinzicrl, Eur. Phys. J. C 41 (2005) 361 
[arXiv:hep-ph/0502165]. 

[15] T. Binoth, J. P. Guillet, G. Hcinrich, E. Pilon and T. Reiter, arXiv:0810.0992 [hep-ph]. 

[16] A. Dcnncr, S. Dittmaier, M. Roth and L. H. Wieders, Phys. Lett. B 612, 223 (2005) 
[arXiv:hep-ph/0502063]. 

[17] R. K. Ellis, W. T. Giele and G. Zanderighi, JHEP 0605 (2006) 027 [arXiv:hep-ph/0602185]. 

[18] A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, JHEP 0808 (2008) 108 
[arXiv:0807.1248 [hep-ph]]. 

[19] W. Giele, E. W. N. Glover and G. Zanderighi, Nucl. Phys. Proc. Suppl. 135 (2004) 275 
[arXiv:hep-ph/0407016]. 

[20] Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Nucl. Phys. B 425, 217 (1994) 
[arXiv:hep-ph/9403226]. 

[21] R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 725, 275 (2005) [arXiv:hcp-th/0412103]. 

[22] R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 715, 499 (2005) [arXiv:hcp-th/0412308]. 

[23] Z. Bern, L. J. Dixon and D. A. Kosowcr, Phys. Rev. D 73 (2006) 065013 
[arXiv:hcp-ph/0507005]. 

[24] C. F. Bcrgcr, Z. Bern, L. J. Dixon, D. Fordc and D. A. Kosowcr, Phys. Rev. D 74 (2006) 
036009 [arXiv:hep-ph/0604195]. 



- 25 - 



C. F. Berger, Z. Bern, L. J. Dixon, D. Forde and D. A. Kosower, Phys. Rev. D 75 (2007) 
016006 [arXiv:hep-ph/0607014]. 

G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. Phys. B 763, 147 (2007) 
[arXiv:hcp-ph/0609007]. 

R. K. Ellis, W. T. Giele and Z. Kunszt, JHEP 0803 (2008) 003 [arXiv:0708.2398 [hcp-ph]]. 

G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0803 (2008) 042 [arXiv:071 1.3596 
[hcp-ph]]. 

C. F. Bcrgcr et al, Phys. Rev. D 78 (2008) 036003 [arXiv:0803.4180 [hcp-ph]]. 

W. T. Gielc, Z. Kunszt and K. Mclnikov, JHEP 0804, 049 (2008) [arXiv:0801.2237 [hcp-ph]]. 

W. T. Giele and G. Zandcrighi, arXiv:0805.2152 [hep-ph]. 

R. K. Ellis, W. T. Gielc, Z. Kunszt and K. Mclnikov, arXiv:0806.3467 [hcp-ph]. 

F. A. Bercnds and W. Giele, Nucl. Phys. B 294 (1987) 700. 

M. L. Mangano, S. J. Parke and Z. Xu, Nucl. Phys. B 298 (1988) 653. 

F.A. Bercnds and W. T. Giele, Nucl. Phys. B 306 (1988) 759. 

Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 437, 259 (1995) 
[arXiv:hcp-ph/9409393]. 

Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 513, 3 (1998) [arXiv:hcp-ph/9708239]. 

Z. Bern, L. Dixon and D. Kosower, Ann. Rev. Nucl. Part. Sci. 46, 109 (1996). 

J. Collins, Renormalization, Cambridge University Press, 1984. 

V. Del Duca, L. J. Dixon and F. Maltoni, Nucl. Phys. B 571, 51 (2000) 
[arXiv:hep-ph/9910563]. 

Z. Bern and D. A. Kosower, Nucl. Phys. B 379 (1992) 451. 

Z. Bern, A. De Freita, L. J. Dixon and H. L. Wong, Phys. Rev. D 66 (2002) 085002. 

David H. Bailey, A Fortran-90 Based Multiprecision System, ACM Transactions on 
Mathematical Software, vol. 21, no. 4, 1995, pg. 379387. 

R. K. Ellis and G. Zandcrighi, JHEP 0802, 002 (2008) [arXiv:0712.1851 [hcp-ph]]. 

J. M. Campbell, E. W. N. Glover and D. J. Miller, Phys. Lett. B 409 (1997) 503 
[arXiv:hcp-ph/9706297]. 

C. F. Bcrgcr et al, arXiv:0808.0941 [hcp-ph]. 

R. Kleiss, W. J. Stirling and S. D. Ellis, Comput. Phys. Commun. 40, 359 (1986). 



- 26 - 



